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Abstract 

We present results from a study of light meson spectra and structure obtained 
within the framework of light-front QCD formulated on a transverse lattice. 
We discuss how imposing Lorentz covariance conditions on meson dispersion 
relations allows determination of parameters in the transverse lattice Hamil- 
tonian. The pion distribution amplitude obtained in this framework is rather 
close to its asymptotic shape. 

I. INTRODUCTION 

Many high energy scattering experiments probe hadron structure very close to the light- 
cone. Hence correlation functions probed in such experiments have a particularly simple 
physical interpretation in the light-front (LF) framework (or infinite momentum frame), i.e. 
in a Hamiltonian framework where x + = x° + x 3 is 'time'. For example, parton distribution 
functions measured in deep-inelastic scattering experiments provide information about the 
LF momentum distribution of quarks in the target hadron. Off-forward parton distributions, 
probed in deeply virtual Compton scattering, have the interpretation of matrix elements of 
the LF momentum density operator between states of unequal momenta. Furthermore, the 
pion distribution amplitude, relevant for example for the asymptotic pion form factor, can 
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be related to the LF wavefunction for the pion in the qq Fock component when q and q have 
the same _L position. 

Even though these important observables have their most physical interpretation on the 
LF, it is, at least in principle, possible to calculate them in any approach to QCD. However, 
what distinguishes the LF framework from all other formulation of QCD is that the above 
observables have a very direct and physical connection to the microscopic degrees of freedom 
in terms of which the Hamiltonian is constructed. Because of this unique feature, it should 
be much easier in this framework to gain a physical understanding between experiment and 
phenomenology on the one hand and the underlying QCD dynamics on the other. 

Other motivations to study QCD formulated on the LF derive from the fact that the 
complexity of the vacuum seemingly shifts from the states to the operators in this framework 
IJ. This results in a separation between the physics of the vacuum and the parton structure 
of hadrons which implies for example that a constituent picture of hadrons has a chance to 
make sense [0. 

Of course, just like in any other approach to QCD, it is necessary to regularize both 
UV and IR divergences before one can even attempt to perform nonperturbative calcula- 
tions. The transverse lattice || is an attempt to combine advantages of the LF and lattice 
formulations of QCD. In this approach to LF-QCD the time and one space direction (say 
x 3 ) are kept continuous, while the two 'transverse' directions xj_ = (x^x 2 ) are discretized. 
Keeping the time and x 3 directions continuous has the advantage of preserving manifest 
boost invariance for boosts in the x 3 direction. Furthermore, since x = x° ± x 3 also re- 
main continuous, this formulation still allows a canonical LF Hamiltonian approach. On 
the other hand, working on a position space lattice in the transverse direction allows one to 
introduce a gauge invariant cutoff on _L momenta — in a manner that is similar to Euclidean 
or Hamiltonian lattice gauge theory. 

In summary, the LF formulation has the advantage of utilizing degrees of freedom that 
are very physical since many high-energy scattering observables (such as deep-inelastic scat- 
tering cross sections) have very simple and intuitive interpretations as equal LF-time (x + ) 
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correlation functions. Using a gauge invariant (position space-) lattice cutoff in the _L di- 
rection within the LF framework has the advantage of being able to avoid the notorious 
1/A; + divergences from the gauge field in LF-gauge which plague many other Hamiltonian 
LF approaches to QCD [§. 

The hybrid treatment (continuous versus discrete) of the longitudinal/transverse direc- 
tions implies an analogous hybrid treatment of the longitudinal versus transverse gauge field: 
the longitudinal gauge field degrees of freedom are the non-compact A^ 1 while the transverse 
gauge degrees of freedom are compact link-fields. Each of these degrees of freedom depend 
on two continuous (x ± ) and two discrete (nj_) space-time variables, i.e. from a formal point 
of view the canonical transverse lattice formulation is equivalent to a large number of cou- 
pled 1 + 1 dimensional gauge theories (the longitudinal gauge fields at each nj_) coupled to 
nonlinear a model degrees of freedom (the link fields) ||. 

For a variety of reasons it is advantageous to work with transverse gauge degrees of 
freedom that are general matrix fields rather than U G SU(Nc)- First of all, we would like 
to work at a cutoff scale which is small (in momentum scale) since only then do we have 
a chance to find low lying hadrons that are simple (i.e. contain only few constituents). If 
one wants to work on a very coarse lattice, it is useful to introduce smeared or averaged 
degrees of freedom. Upon averaging over neighboring 'chains' of SU(Nc) fields one obtains 
degrees of freedom which still transform in the same way as the original SU(Nc) degrees of 
freedom under gauge transformations but are general matrix degrees of freedom no longer 
obeying WU = 1 and det(U) = 1. The price that one has to pay for introducing these 
smeared degrees of freedom are more complicated interactions. The advantage is that low 
lying hadrons can be described in a Fock expansion (this has been confirmed by calculations 
of the static quark-antiquark potential || and glueball spectra [0]). 

Another important advantage of this 'color-dielectric' approach is that it is much easier 
to construct a Fock expansion of states out of general linear matrix fields than out of fields 
that are subject to non-linear SU(Nc) constraints. 

In the color-dielectric approach the complexity is shifted from the states to the Hamil- 
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tonian: In principle, there exists an exact prescription for the transformation from one set 
of degrees of freedom (here U's) to blocked degrees of freedom M = J2av Hi Ui 



The problem with this prescription is that S'e//. is not only very difficult to determine directly, 
but in general also contains arbitrarily complicated interactions. 

A much more practical approach towards determining the effective interaction among 
the link fields nonperturbatively is the use of Lorentz invariance. This strategy has been 
used in a systematic study of glueball masses in Ref. 0, where more details can be found 
regarding the effective interaction. One starts by making the most general ansatz for the 
effective interaction which is invariant under those symmetries of QCD that are not broken 
by the _L lattice. This still leaves an infinite number of possible terms and for practical 
reasons, only terms up to fourth order in the fields and only local (in the _L direction) terms 
have been included in the Ref. ]7|. The coefficients of the remaining terms are then fitted to 
maximize Lorentz covariance for physical observables, such as the QQ potential (rotational 
invariance!) and covariance of the glueball dispersion relation. It should be emphasized 
that these are first principle calculations in the sense that the only phenomenological input 
parameter is the overall mass scale (which can for example be taken to be the lowest glueball 
mass or the string tension) . The only other input that is used is the requirement of Lorentz 
invariance. 

The numerical results from Refs. 0,^] within this approach are very encouraging: 

• with only a few parameters, approximate Lorentz invariance could be achieved for a 
relatively large number of glueball dispersion relations simultaneously as well as for 
the QQ potential 

• the glueball spectrum that was obtained numerically on the _L lattice for Nc — > oo is 
consistent with Euclidean Monte Carlo lattice gauge theory calculations performed at 
finite Nq and extrapolated to Nc — > oo . 




(1.1) 
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For further details on these very interesting results, the reader is referred to Refs. |7||| and 
references therein. In this paper, we are applying this program to mesons on the transverse 
lattice within the femtoworm approximation |J. This Hamiltonian has already been ana- 
lyzed in Ref. |J for mesons with k± — 0. In this work, as well as in a related independent 
paper the work from Ref. || is generalized to mesons with nonzero transverse momenta. 
This generalization allows one to address the issue of Lorentz invariance of meson spectra. 
There are significant differences in the general approach between our work and Ref. ||. 
First of all, while we employ continuous basis function techniques to solve the longitudi- 
nal dynamics, Refs. PJTTj use discrete light-cone quantization (DLCQ). Secondly, there is 
also a difference in the procedure: both this work and Ref. |Tl| use a longitudinal gauge 



coupling taken from pure glue calculations, while Ref. refits this parameter. The latter 
difference has important consequences which we discuss in detail in section III. Our results 



confirm some of the results presented in Ref. [PHI , which serves as a non-trivial test for 
both approaches. However, the main motivation for our work lies in a different direction. 
Our goal was to understand to what extent the requirement of imposing Lorentz invariance 
sufficiently constrains the parameters of the Hamiltonian first studied in Ref. [§] and also 
to investigate whether this Hamiltonian and the femtoworm approximation are able to yield 
Lorentz invariant results. The paper is organized as follows. We briefly describe the Hamil- 
tonian, the approximations used to solve it and some of it's subtle properties in section II. 
In section III, we provide a perturbative analysis of the tt — p splitting. The results from 
this section will be of importance for understanding the non-perturbative numerical results 
presented in section IV. This is followed by some comments and comparisons of this work 
with another concurrent work in section V and a brief summary in section VI. 

II. THE HAMILTONIAN 

The transverse lattice Hamiltonian for Wilson fermions that forms the basis of our work 
and Ref. [JTTJ] was first constructed in Ref. [§], and we refer the reader to these works for 
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more details. Here we restrict ourselves to a more qualitative description of the degrees 
of freedom that enter the Hamiltonian and we follow up later with more specifics when 
we discuss the integral equations describing mesons in the femtoworm approximation. The 
degrees of freedom that enter the Hamiltonian are quark and antiquark creation operators 
at each site and link-fields (which, in the spirit of the color-dielectric approach ]7|], we take 
as general matrix fields) on the transverse links connecting the sites. For simplicity, we 
work in the large Nc limit. The first approximation that we use is a light front Tamm- 
Dancoff truncation of the Fock state expansion upto and including the 3-particle sector. 
Upon reminding oneself that the fermion degrees of freedom (quarks and antiquarks) always 
occupy the lattice (transverse) sites and the gauge degrees of freedom (link fields) reside on 
the lattice spacings, it is easy to understand the following two direct consequences of the 
aforementioned truncation: 

• The 2-particle states consist of a quark and an antiquark sitting on the same lattice 
site. 

• The 3-particle states consist of a quark and an antiquark separated by at most one 
link field. 

Since link fields always reside on the lattice spacings as demanded by gauge invariance, the 
3-particle Fock space is therefore spanned by all different configurations of a quark and an 
antiquark separated by one link field. This dumb-bell shaped configuration of the 3-sector 
naturally requires one to introduce an index to specify it's transverse orientation. We do so 
as shown in Fig. |l[ 

In keeping with the truncated Fock space, one can therefore write the meson state vector 

as: 

| meson) = c<i \ qq) + C3 | qqg) (2.1) 

The total angular momentum of the fermions in the 2-sector is contributed wholly by their 
respective spins since there is no contribution from orbital angular momentum when they 
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are both on the same site. On the other hand, the total angular momentum of the 3-particle 
configuration is the sum of both their spins as well as the orbital angular momentum typical 
of dumb-bell shaped systems. To characterize a 2-particle state, one needs the spins of the 
two quarks along with their momenta, both longitudinal as well as transverse. Additionally, 
a 3-particle state will require the direction of it's configuration (see Fig. [TJ) in the discretized 
xy directions. The Hamiltonian, however, remains invariant when translated in multiples 
of the transverse lattice spacing. Consequently, the transverse momenta of the mesons 
are conserved (modulo it). This allows one to solve the light front Hamiltonian equation, 
Hlf I ip) — P I for fixed values of the transverse momentum, k±. For a given value of 
k± we therefore make the following ansatz: 

ft± 

dl(l-x,n ± )i> sl>S2 (x) |0) (2.2) 



qqg) = E e 



1 rl-x 



dx dybl i (x,n ± )dl 2 {y,n ± + 



o Jo 

a f (l - x - y)ip SliS2 (x, y, e h ) | 0) (2.3) 

where the spin degrees of freedom, s\, S2 E (Til) an d the orbital degrees of freedom, e*j E 

(±e x ,±e y ). In the above, ^ SljS2 (x) and ip Sl>S2 (x ) y ) e i ) are the wavefunction amplitudes in 

the 2-particle and 3-particle sectors respectively. It is a straight forward though tedious 

procedure to show that the QCD Hamiltonian with the above state vectors in the large N c 

limit and under the light front gauge yields the following integral equations in terms of the 

wavefunction amplitudes in the 2-particle and 3-particle sectors respectively: 

(rrA — l m 2 g - 1\ 
y x 1 — x J 

1 dx' ^ 81 ' 82 ^ (2.4) 
(x — x') 2 



fax,* (x, y, et) - / dx — — === 

J ° \x - x'Y 2^x u x' u 

l ~ x d ,Tj>ai,a a (x,y r ,g i ) yu + y'u 







(y - v'Y 2^yl 

(2.5) 

In Eq. ( p.4|) , the first two terms on the right hand side are the kinetic energies of the quark 
and the anti-quark while the third term describes the Coulomb interaction between them. 
In Eq. (|2.5| ), the first three terms on the right hand side are the kinetic energies of the 
quark, anti-quark and the link fields, the fourth term describes the Coulomb interaction 
between the quark(antiquark) and the link field and the fifth term describes the Coulomb 
interaction between the antiquark(quark) and the link field. It should be noted that a 
Coulomb interaction term between the quark and antiquark in the 3-particle sector does not 
appear since it is disallowed under the large N c limit. Three remarks are in order at this 
point. First, if 22 and if 33, i.e. the terms in the Hamiltonian that are diagonal in Fock space 
are independent of the quark spin and the orientation of the link. Second, the ^'s in Eqs. 
( j2.4| ) and fl2.5p still depend implicitly on k±_ because although k±_ is conserved, k±_ still enters 
the integral equations as a parameter. Finally, Eqs. ( |2.4|) and (|2.5| ) have been rescaled in a 
manner such that G 2 = 1. 

The natural question that arises at this point is: How do the 2-sector and the 3-sector 
mix with each other. This brings us to the issue of propagation of the mesons on the lattice. 
Qualitatively, the way this takes place is as follows: A q(q) hops to a neighboring site with the 
emission of a link field (2-sector — > 3-sector). The q(q) sitting on the original site then jumps 
to join the q(q) by absorbing the link field (3-sector — > 2-sector). Since this propagation 
resembles that of an inchworm on a scale of a femtometer, we were inspired to called this the 
femtoworm approximation. The hopping of a quark or an antiquark from one site to the next 
can cause it to either flip its spin or not. Accordingly, we introduce two distinct couplings 
between the fermion and the link degrees of freedom. In this work, we label the spin-flip 
coupling as m v and the non spin-flip coupling as m r . The corresponding interaction terms 
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are presented in details later in this section. The natural question that now arises from the 
above is: Is it possible for a link field that is emitted by a quark(antiquark) through the 
m r coupling to be absorbed by an antiquark(quark) through the m v term, in other words, 
'mixed' hopping ? The answer to this is 'No'. To understand the hopping interactions better 
one needs to analyze the rotational properties of the Hamiltonian more closely. For = 0, 
the _L lattice Hamiltonian is invariant under rotations around the z — axis by multiples of 
7r/2, giving rise to conservation of total angular momentum modulo 4. Since J z = S z in the 
2 particle sector and S z can only assume the values —1,0, and 1, S z in the 2 particle sector 
is conserved. Since mixed hopping would change S z by one unit, this means that the sum 
of all contributions from mixed hopping must add up to zero. The observation that helicity 
conservation (modulo 4) for k± = prohibits interference between the two hopping terms 
was already made in Ref. ||. 

The interaction of the 2-particle sector with the 3-particle sector will therefore have four 
different possibilities: 

• a quark(or an antiquark) in the 2-sector emits a link through the m v term and flips 
it's spin. 

• a quark(or an antiquark) in the 2-sector emits a link through the m r term and does 
not flip it's spin. 

Of course, the interaction to go from the 3-sector to the 2-sector has the opposite effect, i.e, 

• a quark(or an antiquark) in the 3-sector absorbs a link through the m v term and flips 
it's spin. 

• a quark(or an antiquark) in the 3-sector absorbs a link through the m r term and does 
not flip it's spin. 

In order to understand the spin dynamics due to the hopping terms as described above we 
present below the interaction Hamiltonian in all its glory. 
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Spin flip hopping in the x-direction: 

8Hl Up = 




-b\(n ± ,pj)b^(n ± + e x ,pf) 
-b\(n ± ,p~j)b l (n ± - e x ,pf) 
+b\(n ± ,p})bi(n± - e x ,pf) 
-d{(n ± ,p^)d|(n ± + e x ,pf) 
+d\(n ± ,pj)d^(n ± + e x ,pf) 
+d\(n ± ,pj)d l (n ± - e x ,pj) 
-d\(n ± ,pj)d^(n ± - e x ,pf) } 

Spin flip hopping in the ^/-direction: 

5H f y lip = 

\Pf Pi J n± 

+b\(n ± ,pf)bi(n± + e y ,pf) 
-fc{(ri±,p+)^(n± - e y ,p+) 
-b\(n±,pf)bi(n± - e y ,p+) 
-d\(n ± ,ppd l (n ± + e y ,pf) 
-d\(n ± ,ppd 1 (n ± + e y ,pj) 
+d\(n±,pj)d l (n± - e y ,pf) 
+d\(n ± ,pj)d^(n ± - e y ,pf) } 

Hopping without spin flip: 

- 5 + 5 )T,[ b }( fl ^P^)h(n± + e x ,pf) 

\Pf Pi ) nj 



+b\(n ± ,ppb r (n± - e x ,pf) 

+b\(n±,pj)b l (n± - e x ,p+) 
+d\(n±,p~})d< [ (n± + e x ,pf) 
+d\{n±,p})d<t(n± - e x ,pf) 
+d\(fi±,pj)d i (n± + e X} pf) 
+d\(n±,pj)di(n± - e x ,pf) 
+b\(n ± ,p^)bi;(n ± + e y ,pf) 
+b\(n ± ,p^)bi;(n ± - e y ,pf) 
+b\(n ± ,p~})b l (n± + e y ,pf) 
+b\(n ± ,p})b l (n ± - e y ,pf) 
+d\{n±,ppdf(n± + e y ,pf) 
+d\{n±,Pf)df(n± - e y ,pf) 
+d\(fi±,pj)d i (n± + e y ,pf) 
+d\(n ± ,p^)d l (n ± - e y ,pf) } 

Implicit in the above expressions are link field creation and annihilation operators on the 
link connecting the initial and final state. 

The transition matrix elements are obtained by sandwiching the above interactions be- 
tween the 2-particle and the 3-particle states as written in Eqs. (|2.2|) and (|2.3|) and are 
provided in the table in appendix B. 

The schematic structure of the Hamiltonian therefore looks like the following: 



H 







H22 


H23 




ip(x) 






H32 


H33 




4>(x,y) 



(2.9) 

where if 23 and H 32 involve all the transition matrix elements. 
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In order to evaluate all the matrix elements, we expanded the wavefunction amplitudes 
using a basis of continuous functions as follows: 

A 



k=0 1=0 



(2.10) 



(2.11) 



where the A in the Eq. 2.10 and that in Eq. 2.11 need not be equal. 

In the table (see appendix B), F r q (F|) are the overlap integrals when a quark(antiquark) 
in the 2-sector goes to the 3-sector by emitting a link field through the m r hopping term 
(the spin non-flip term) and f q (fq) are the overlap integrals when a quark(antiquark) in 
the 2-sector goes to the 3-sector by emitting a link field through the m v hopping term (the 
spin-flip term). See appendix A for more details of these overlap integrals. 

A factor (jpj was introduced in the integrand of the transition matrix elements to 
regularize the infrared (small x, x') divergence that arises from the self energy of quarks 
in the case when the quark that emits a link field also absorbs it. Of course, all physical 
observables will have to be independent of the value of e as will be demonstrated in section 
IV. 

The result of all the above is the following matrix equation for each k and k, I: 



H 



N k , k > 
N k ^ k > v 



a k > 
a k 'v 

M kk i Mk,k'i' 

Mk'i>,k Mki,k'i' 



a k > 
a k 'i' 



(2.12) 



where N ktk i and Nki^'v are the norm matrices, M kk i and Mki t k'i' are the matrix elements in 
the pure 2-particle and 3-particle sectors, respectively and finally M^k'i 1 an d Mk>i',k are the 
transition matrix elements. The above matrix equation is then converted into an eigenvalue 
equation that yields eigenvectors in terms of the a k s and a^s. 



12 



Another significant result of this approximation that is worth mentioning at this point 
is that unlike the low lying mesons, the higher excited states are not expected to show 
the correct dispersion. The reason for this is that since the higher states are very excited, 
the quark and the antiquark will be physically farther apart. This would mean that these 
mesons will mostly or always remain in the qqU (where U is the link field) state (see Fig. 
H) and for it to propagate, it will need to go to a qqUU state or states with larger number 
of link fields. But since our approximation does not allow states with more than one link 
field, these excited states simply cannot propagate. 

III. PERTURB ATIVE ANALYSIS OF THE ir-p SPLITTING 

In order to gain a qualitative understanding about the interplay between different pa- 
rameters in the _L lattice Hamiltonian, it is instructive to study a simple model, where one 
treats the admixture of the 3-particle Fock component to the tt and p as a perturbation. 

To th order, i.e. when the coupling between 2 and 3 particle Fock component is turned 
off, there is no spin dependence of the interactions and the tt and p are degenerate. Likewise, 
there is no 'hopping' (i.e. _L propagation) of mesons and thus energies are independent of 
kj_ giving rise to an infinite _L lattice spacing (in physical units). 

In the next order we treat the coupling between 2 and 3 particle Fock components as a 
perturbation (note that interactions which are diagonal in the particle number, such as the 
confining interaction in the longitudinal direction are still treated non-perturbatively). As 
described in section II, there are two interactions that mix Fock sectors: hopping due to the 
m r -term (without helicity flip) and hopping due to the vector coupling m v (with helicity 
flip). Starting from a basis of ! t Hooft eigenstates which are plane waves in the _L direction 
and where the qq in the 2 particle Fock component carry spins | jf), | ][), | ||), and | [[) 
respectively, one thus finds for the energy in second order perturbation theory to be: 
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H=M$- Ml 



C-x ~t~ Cy 
C x + C 










V 



+ M 2 . 











Cy c x 









H~ Cy 











H~ Cy 

Cx H~ 





Cx Cy 






Cy C x 









(3.1) 



Here M\ r and M 2 ^ are some second order perturbation theory expressions involving 
matrix elements between 2 and 3 particle states that are eigenstates of the diagonal parts 
of the Hamiltonian (kinetic + Coulomb), and c« = cosfcj. 

Several general and important features can be read off from this result. First of all, and 
most importantly, Eq. ( |3.1| ) shows that the r-term gives rise to a dispersion relation with 
the same _L speed of light for the 7r and the p's, while the vector interaction breaks that 
symmetry. This observation already indicates that it may be desirable to keep the r-term 
much larger than the spin-flip term. We will elaborate on this point below. 

At k x = k y = 0, the eigenstates of the above Hamiltonian are the p±i i.e. | ft) and 



II), with M 2 = M. 



±i 



Ml 



M\ r , the |p ) = | tl + It), with M 2 = M± x + M\ v and 



the 1 7r) = | tl — IT), with M 2 = M± x — M 2 V For nonzero _L momenta, there will in general 
be mixing among the p + i and the p_i, but not among the other states since helicity in the 
2-particle Fock sector is still conserved modulo 2. Expanding around = 0, and denoting 
M 2 = Mq — M\ r one finds to Oik'jJ the following eigenstates and eigenvalues. 
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state 


IVl \y J 


M 2 (\r 2 

IVl 


) - M l u J 


1 1 it 


M - M l v 


(.2, 1,2 
T\/r2 K x+ K y 

M l,r 2 


1.2 ,1,2 
+ M l,v 2 


ti i it 
H + ll 


i rl i i /2 
M + M lv 


M l,r 2 


M l,v 2 


TT-II 


M 2 


„ 1,2, u2 
1\/T2 K x + K y 

m l,r 2 


71 /T2 k x K y 
1V1 l,v 2 


TT + U 


M 2 


A/T2 K x+ K y 

IVI l,r 2 


1 /i #2 K y 
+ 2 



(3.2) 



Eq. ( ^.2[ ) illustrates a fundamental dilemma that hampers any attempt to fully restore 
Lorentz invariance within the femtoworm approximation: M\ v not only governs the splitting 
between the 7r and the (h=0) p but is also responsible for violations of Lorentz invariance 
among the different helicity states: If one determines the _L lattice spacing in physical units 
for each meson separately, by demanding that the _L speed of light equals 1, one finds for 
example 



1 
1 



Ml + Ml 
Ml - Ml 



(3.3) 



i.e. increasing the 7r— p splitting is typically accompanied by an increase in Lorentz invariance 
violation 



1 



M 2 — M 2 . 



(3.4) 



For the p±\ the breaking is of a similar scale, plus one also observes an anisotropy in the 
dispersion relation on the same scale. 

Therefore, in order to avoid a large breaking of Lorentz invariance, it will be necessary 
that 



Ml . > M 2 o - Ml 



(3.5) 
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If one keeps the it — p splitting fixed at its physical value then there are two ways to achieve 
this condition. One possibility is to simply increase the Yukawa coupling that appears in 
the r-term. This increase of the r-term tends to decrease the _L lattice spacing for both 7r 
and p's (see Fig. 0) and in order to achieve satisfactory Lorentz invariance (in the sense of 
uniform _L lattice spacings) one needs to make the lattice spacing smaller than the Compton 
wavelength of the p meson. However, one cannot make the r-term coupling arbitrarily large 
because at some point there occurs an instability (tachyonic M 2 !). Such instabilities for large 
coupling are common in the LF formulation of models with Yukawa coupling and might be 
related to a phase transition (similar to the phase transition in <j) 4 theory that occurs as the 
coupling is increased). 

Fortunately, there exists another possibility to make these matrix elements large, without 
increasing the Yukawa couplings. This derives from the fact that the hopping interactions are 
proportional to ( \ ± -^p, where x (x') are the momenta of the active quark before(after) 
the hopping. Because of the singularity as x, x' —>■ 0, matrix elements of the hopping terms 
are greatly enhanced if the unperturbed wave functions are large near x = and x = 1. 
Since the unperturbed wave functions in the 2 particle Fock component vanish like x 13 near 
x = 0, where (5 oc matrix elements of the hopping interaction become very large when 
one makes ^ very large. 

Therefore, the larger one chooses the more one restores Lorentz invariance of the it 
and p dispersion relations because one can keep the ir-p splitting fixed while decreasing the 
coupling of the spin flip interaction. At the same time, keeping the r-term interaction fixed 
one increases the dominance of the r-term contribution in \ and thus not only reduces the 
lattice spacing in physical units, but also obtains dispersion relations for the n and the p's 
that look more and more similar — as demanded by Lorentz invariance. Unfortunately, we 
are not completely free to pick whatever value of ^ we like because m 2 q and G 2 are largely 
fixed by the center of mass in the 7r-p system as well as by fitting the physical string tension 
in the pure glue sector. 

The above analysis may also shed some light on the shape of the pion distribution 
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amplitude in Ref. ||, which is much more flat (actually, it even exhibits a slight double- 
hump shape) compared to our results (see next section). The reason is that the numerical 
value for G 2 in Ref. |J was not taken from the pure glue calculation, but was rather used as 
a free parameter, which was also varied in order to maximize Lorentz covariance. Although 
such a first principle approach would normally be desirable, minimizing Lorentz covariance 
violations for the ir and p dispersion relations drives G 2 — > oo. When the longitudinal gauge 
coupling G 2 is very large compared to other scales in the problem, the wavefunction tends 
to become very flat — a phenomenon familiar from the 't Hooft equation. Of course, since 
a DLCQ cutoff was employed in Ref. |J, G 2 never really become infinite, but through the 
extrapolation procedure the near flat shape for the pion distribution amplitude emerged. In 
a more recent paper [|TTJ, where Dalley adopted a procedure similar to ours and G 2 was kept 
fixed, a pion distribution amplitude is obtained that is much more consistent with ours. 



IV. NONPERTURBATIVE NUMERICAL RESULTS 

One cannot expect this Hamiltonian to describe the higher excited states because of the 
truncation of the Fock space upto the 3-particle sector. This is demonstrated in Fig. |3| 
where we show typical results for the probability to find each state to be in the 3-particle 
sector. One finds that whereas the low-lying states have a fairly large probability to be in 
the 2-particle sector, it is quite the opposite for the higher states. We therefore expect that 
excited states are much more strongly affected by the ommission of 4 and more particle states 
and consider in the following only the 7r and p states in our study of Lorentz covariance. 

The parameters in the Hamiltonian are: 

• longitudinal gauge coupling 

• r-term coupling (hopping without spin flip) 

• spin-flip coupling (hopping with spin flip) 

• kinetic mass (2 particle sector) 
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• kinetic mass (3 particle sector). The observables that we studied showed little depen- 
dence on this parameter, so we kept it fixed at a constituent mass. 

• link field mass. Similar to kinetic mass in 3 particle sector. Furthermore, demanding 
Lorentz invariance in the pure glue sector, one finds a renormalized trajectory. This 
makes sense since the _L lattice scale is unphysical. We keep the link field mass fixed 
at a value which yields relatively small _L lattice spacing. 

In the spirit of the color-dielectric approach J7|, we would like to use only chiral symmetry 
(small pion mass), Lorentz invariance and one length scale (e.g. physical p mass or string 
tension) as the only input into our calculations. However, as explained in the previous 
section, if one uses as input only the n and p masses, demanding exact Lorentz invariance 
within the framework of the femtoworm approximation would drive G 2 /m 2 p to infinity and 
would thus give rise to string tensions that are inconsistent with phenomenology. Likewise, 
as input only the ir mass together with Lorentz invariance would give rise to a degenerate 
7t-p system. It is thus clear that, as long as one does not go beyond a truncation of the Fock 
space above 3 particles (q, q and one link quantum), a first principle calculation (where only 
one mass scale is used as input) will not be possible. In the following, we will therefore use 
two phenomenological scales as input parameters: the physical value of the p mass, as well 
as the physical string tension. 

Although the dependence of physical parameters on the input parameters is in general 
rather complex and non-perturbative, one can understand at a qualitative level how the 
input parameters influence the relevant physical scales: The physical string tension in the 
longitudinal direction determines G 2 in physical units. 

Modulo mass renormalization due to the Yukawa couplings, the quark mass in the 2- 
particle Fock component and the gauge coupling are strongly constrained by fitting the 
physical string tension (which determines G 2 and the center of mass of the n-p system. This 
leaves us with only the r-term and the helicity flip hopping term couplings as parameters to 
vary. The helicity flip term is not only responsible for n-p splitting within our approximations 
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but also for violations of Lorentz invariance (different _L lattice spacings in physical units 
for different mesons). 

The observables that we studied, including Lorentz invariance, show rather little sensi- 
tivity [] to the precise values of the 3-particle sector quark masses (Fig. |]), which we thus 
keep at a value corresponding to a constituent mass (about half the p mass). This leaves us 
with 4 free parameters. 

Using a value of G 2 = 9 ^° ~ 0.4 GeV 2 in physical units, as determined from the 
string tension in Ref. (which is larger than the one previously used [§,0]), we were able 
to produce the physical n-p splitting with only a relatively small spin-flip coupling. This 
allows us to chose the r-term large enough so that r-term hopping dominates over spin-flip 
hopping and therefore violations of Lorentz invariance (as measured by comparing quadratic 
terms in the dispersion relation of mesons in the n-p sector) are only on the order of 20%. 
The average _L lattice spacing (from the dependence of the energy on Pj_) is found to be 
aj_ w 0.5/m. For the 7r wave function, we find a shape (see Fig. |5]) that is very close to the 
asymptotic shape 4> asy (x) = 6x(l — x). 

This is surprising if one considers that the lattice spacing is still relatively large and 
hence the momentum scale is still very low. The p meson distribution functions are shown 
in Figs. § and |7|. Notice that the momentum dependence of the ir dispersion relation is 
stronger than the one for the p and for the parameters used in our calculations this even 
leads to a level crossing f] The reason for this behavior is very simple. As pointed out in 
Ref. ||, for k± = spin-flip hopping is attractive for the 7r and repulsive for the po, while 
the opposite is true at the end of the Brillouin zone. More precisely, for k x = 0, the hopping 
term in the x direction is attractive for the 7r and repulsive for the p° with the opposite 
being true for k x = 7T. An analogous statement holds for the ^/-direction. If there were no 



As long as we kept the other parameters floating! 

2 Note that there is no mixing at the crossover point since the tt and p have opposite C-parity. 
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r-term present, the dispersion would arise solely from the spin flip term and the po in the 
corner of the Brillouin zone would be degenerate with the tt at k± = — a phenomenon 
which is nothing but a manifestation of species doubling at the hadronic level. Of course, 
since we have included an r-term, the p° in the corner of the Brillouin zone is no longer 
degenerate with the tx at k± = 0, but the level crossing remains (see Fig. ^). 

Similarly one can understand the degeneracy of the 7r and p° at the edge of the Brillouin 
zone k± = (tt, 0). At that point, the x-hopping term is repulsive for the tt and attractive for 
the p . The opposite is true for the y- hopping term and the two contributions cancel one 
another leaving the tt and the p° degenerate. 

From the normalization of the tt distribution function we were able to extract a numerical 
estimate for the tt decay constant f^p^ = (OlA^lp) in terms of the H — helicity component 
of the tt wavefunction in the two particle Fock component $ + _(x) 

2 [Nr r 1 

U = —\ — / dxip + ^(x) « 300 MeV, (4.1) 
a± V tt Jo 

which is much larger than the experimental value f* xp ~ 93 MeV. We expect that including 
higher Fock components will reduce the size of the two particle Fock component and hence 
also f w . 

The parton distribution function for the tt is shown in Fig. The fact that q(x) goes 
to zero at small x is an artifact of our e-cutoff and we expect a nonzero value for q(0) in the 
femtoworm approximation as e — ► 0. 

Finally, we also calculated the electro-magnetic form factor for the 7r for purely transverse 
momentum transfers (where overlap formulas can be used) and we find an rms radius for 
the 7i of r\ = —6^2 F(Q 2 )\q =q w (0.243 fm) 2 . Obviously, since our lattice spacing is only 
about 0.5 fm and since we do not let the quark- ant iquark separation to exceed one lattice 
unit, one should not expect very good results for this observable. Since the rms-radius is 
already poor, displaying the form factor is not very meaningful. 

As mentioned in section II, physical observables should be independent of the infrared 
cut-off, e (as long as the cutoff is small). We verified this cut-off independence for a variety 
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of observables. As an example, we demonstrate in Fig. || the cutoff independence of the 
dispersion relations for light mesons. The dashed line is for e = 0.05 and the continuous line 
is for e=0.1. 



V. COMMENTS AND COMPARISONS 

We feel it is worthwhile to compare our results with those in a concurrent paper that 
appeared as we were completing the manuscript for this paper. There are a few noteworthy 
technical differences between the two. Firstly, the treatment of the longitudinal dynamics 



are different (DLCQ in Ref. |TTJ while we employed continuous basis functions). Secondly, 



the pure glue string tension are different (~ I fm vs. ~ 0.5 fm) Otherwise, both this work 



and Ref. [[0] were based on the formalism outlined in Ref. 



We used continuous basis functions because there are some unresolved issues regarding 



the convergence of DLCQ calculations in terms of the DLCQ resolution [|T^]. However, 
because the numerical approximations in both works are very crude, these issues are not 
yet very important and therefore the main difference in this context is more in calculational 
details. 

Because both works used the femtoworm approximation, the rms radius of the 7r is rather 
low compared to the experimental value. The lattice spacing, which defines the maximal 
gg-separation, is slightly larger in Ref. (Tl|. Consequently, the resulting rms radius in Ref. 
|TT| , though still much smaller than the experimental result, is slightly larger than ours. 



Here we should also point out that the algebraic expression used in Ref. [ I 1. 1 for calcu- 



lating f w is incorrect. With our normalization, the result for f w in Ref. \\ l.| increases by a 



factor of three |T5[ and becomes much more similar to ours. 

Since both works determine the parameters in the effective theory by imposing Lorentz 
symmetry as a nonperturbative renormalization condition, we feel that it is important to 
have a handle on the violation of Lorentz invariance (or lack thereof) within the femtoworm 
approximation. That this violation is inevitable was first pointed out in Ref. [H (in which 
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some results from this work were already presented) and has been further reiterated in 
section III of this paper. Since it is impossible to achieve both Lorentz invariance and the 
correct tc — p splitting simultaneously, one, therefore, needs to have a quantitative idea of the 
degree to which Lorentz invariance can be achieved while simlutaneously yielding the correct 
7r — p splitting. Although Ref. |]TT| provides several x 2 plots, it is not quantitatively clear 
from these plots how large the violations of Lorentz invariance are for the it — p dispersion 
relations at the minimum of \ 2 ■ It is therefore not possible for us to compare the two works 
with resepect to the degree to which Lorentz symmetry is violated. Here we should also 
point out an interesting connection between this work and Ref. on the one hand and 
Ref. || on the other: although the spirit of Ref. |J was more phenomenologicalf], the final 
results (parton distribution functions, distribution amplitudes) obtained in Ref. |8j were very 



similar to both Ref. [11] and this work. 



Another difference between our results and those in Ref. Ill is that our 7r distribution 



amplitude is closer to the 'asymptotic' shape than the one in Ref. [jnj - a behavior which 
is also reflected in the parton distribution function. However, this difference is not really 
significant and it may just be due to the differences in the handling of the longitudinal 
dynamics which may favor slightly different sets of parameters. 



VI. SUMMARY AND CONCLUSIONS 

We presented results from a study of light mesons in QCD formulated on a transverse 
lattice using light-front quantization. In transverse lattice QCD the gauge field in the 
(continuous) longitudinal directions is described by a non-compact vector potential, and in 
the (discrete) transverse directions by a non-compact link field. Our approach was based on 



3 The quark mass in the three particle Fock component in Ref. || was simply taken to be a 
constituent mass rather than being fitted to optimize Lorentz symmetry and the ratio of the flip 
to nonflip coupling was simply kept fixed in the fit of mass spectra. 
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the color-dielectric formulation 0, where one assumes that the fields have been "smeared" 
over a finite area and thus are no longer subject to a strict SUc{N) constraint. 

The smeared degrees of freedom are subject to effective interactions which reflect the 
dynamics of the underlying microscopic degrees of freedom. Although this effective interac- 
tion is not known a priori, it is highly constrained because it must still be gauge invariant 
and it should still give rise to Lorentz invariant observables. In this work we investigated 
transverse lattice QCD with fermions. The starting point was an ansatz for the effective 
interactions which contained interaction terms up to fourth order in the fields which are still 
invariant under those symmetries not broken by the transverse lattice. The goal of this work 
was to determine the coefficients in this ansatz by demanding enhanced Lorentz symmetry 
from physical observables. As input parameters we used the physical masses of the 7r and 
the p as well as the string tension. What we found is that the splitting between the 7r and p 
mesons strongly constrains the quark helicity flip interaction term in the LF-Hamiltonian. 
Meson propagation in the _L direction involves helicity flip and nonflip hopping terms for 
the quarks. Once the nonflip hopping term was constrained by fitting the observed it — p 
splitting, the only remaining parameter on which the physical lattice spacing (as determined 
from the meson dispersion relation) depends very strongly is the helicity nonflip term. Even 
though we scanned the parameter space thoroughly, we were only able to find matching val- 
ues for a± for the 7r and p mesons within about a 10% - 20% error. Similar Lorentz invariance 
violations were found in the splitting between p mesons of helicity zero and one. We believe 
that both these errors are due to the truncation of the Fock space (which may contain an 
infinite number of particles) as well as the ansatz for the effective interaction (which may 
contain an infinite number of terms). In fact, given these rather crude approximations, it 
is actually remarkable that the violations of Lorentz invariance in the light meson sector 
are only about 10-20% . The model seems very promising though it displays some basic 
limitations as far as restoration of full Lorentz invariance is concerned. This problem was 
first identified in one of our earlier papers [TIJ but subsequent new understanding (T2j about 



the couplings helped isolate the problem and reduce it significantly for the low-lying mesons. 
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For the higher mesons, the femtoworm approximation is not expected to work unless more 
link field degrees of freedom are introduced, a task that is possible but computationally 
non-trivial. 

While we were completing this work, similar studies of mesons on the transverse lattice 
appeared P,P"T|, which were based on DLCQ rather than basis function techniques for han- 



dling the longitudinal dynamics. However, there is another important difference between 
our work and Ref. ||: we took the numerical value for the longitudinal gauge coupling G 
that was obtained in pure glue studies on the _L lattice [0], while Ref. || allowed G to 'float' 
in addition to the other parameters in order to restore Lorentz symmetry. As we explained 
in Section III, increasing G, as compared to other dimensionful parameters, tends to de- 
crease violations of Lorentz symmetry in the n/p dispersion relations. Therefore, the results 
presented in Ref. @ where extrapolated from calculations performed with values of G that 
were much larger than ours (and larger than G used in the pure glue calculation in Ref. J7j). 
This difference in procedure explains to some extent the difference in Results between this 
work and Refs. ||: most importantly, if one increases G in physical units, while keeping 
the physical pion mass fixed, tends to decrease all other dimensionful parameters The 
pion distribution amplitudes presented in Ref. 0, where G was allowed to float completely 
freely, although slightly 'double humped', is actually numerically consistent with an almost 
completely flat distribution since it was obtained by extrapolation in the DLCQ parameter 
and there is some uncertainty in this procedure, and this is consistent with the above ob- 
servation regarding G. In Ref. [|ll[j, G was no longer allowed to float and the results for the 
distribution amplitude are much more similar to ours. 
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4 One can understand this fact by comparing with the simpler but similar 1+1 dimensional 't 
Hooft model 
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APPENDIX A: OVERLAP INTEGRALS 

To understand the overlap integrals, let us look at two specific elements of the table (see 
appendix B), namely, the first two elements of the first row. For the first element, clearly 
there is no spin flip and the hopping is along the y-direction. There are two possibilities 
through which this transition can take place. Pictorially, they are as shown in Fig. [TTJ. 

Notice that both diagrams yield the same 3-sector lattice configuration which is shown to 
the right. These two diagrams correspond to the terms containing b^(n±,p^)b^ (n± — e y ,pf) 
and rf|(nj_, p~f )d^(n± + e y ,pf) in the the spin non-flip interaction Hamiltonian. When these 
two terms are sandwiched between | qq) and | qqg) from eqns. 2.2 and 2.3, one obtains the 
first matrix element of the table. In particular, and F| are obtained when the integrals 
over the quark and the antiquark momenta are carried out in the above matrix element as 
follows: 



Fl = Jdxjdynx,y)(^ + ^ 



sjl-x-y 
F r q =JdxJ dyr(x,y)( T ^ + - 



yj 

i/>(l-y) (Al) 



y i — x 
1 



ip(x) (A2) 



V 1 - x - y 

Similarly, let us look at the second element of the first row of the table. Clearly, this 
is a spin flip interaction and the hopping takes place in the ^/-direction. Pictorially, this is 
shown in Fig. [1~2| . 

This diagram corresponds to the term in SH^ hp containing c?|(nj_,pj )di(n±+e y ,pf). The 
equivalent 3-sector lattice configuration is once again shown to the right. The corresponding 
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matrix element is therefore obtained by sandwiching this term of the Hamiltonian between 
| qq) and | qqg) and in particular fg is obtained when the integrals over the quark and the 
antiquark momenta are carried out in the above matrix element as follows: 



f q = J dx J dyip*(x, y) I ^ - 



y/l-x-y 



y ) 

V>(i-y) (A3) 



f q = I dx I dyip*(x,y) 



y l — x J 

1 tl){x). (A4) 



y/l-x-y 
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APPENDIX B: TRANSITION MATRIX 
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FIG. 1. Four possible orientations of the 3-particle states. 
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FIG. 2. Dependence of a± on the non-spin flip coupling. 
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FIG. 3. Probability of states to be in the 3-particle sector. 
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FIG. 4. Dependence of the _L lattice spacings on the masses in the 3-particle Fock component 
(other parameters are adjusted correspondingly to keep the physical n and p masses fixed). 
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FIG. 5. Light-cone momentum distribution amplitude for the ir obtained on the _L lattice. The 
dashed line is the asymptotic shape. 
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FIG. 6. Light cone distribution amplitude for p° obtained on the _L lattice. 
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FIG. 9. e independence of dispersion relations. The dashed line is for e = 0.05 and the contin- 
uous line is for e=0.1. 
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FIG. 10. Pion parton distribution function. 
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FIG. 11. Spin non-flip hopping. 



40 



X 



X 




m. 



n± + e % 



V on ± 



FIG. 12. Spin flip hopping. 
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